Coalescence in low-viscosity liquids 
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The expected universal dynamics associated with the initial stage of droplet coalescence are 
difficult to study visually due to the rapid motion of the liquid and the awkward viewing geometry 
Here we employ an electrical method to study the coalescence of two inviscid droplets at early 
times. We measure the growth dynamics of the bridge connecting the two droplets and observe a 
new asymptotic regime inconsistent with previous theoretical predictions. The measurements are 
consistent with a model in which the two liquids coalesce with a slightly deformed interface. 

PACS numbers: 47.55.df, 47.55.D-, 68.03.-g, 47.55.nk, 47.55.N- 



When fluid drops merge, a dramatic transformation 
occurs: the topology changes as the fluid masses, origi- 
nally separated, merge into a single entity. At first, the 
drops are separated by only a small distance. Then a 
thin fluid bridge is formed between them which rapidly 
widens due to surface tension forces. We employ an elec- 
trical method to explore drop coalescence in a low viscos- 
ity fluid and find that, shortly after its initiation, there 
is an unexpected new regime that is dominated by the 
overall deformability of the drops. 

Soon after the instant of coalescence, the length scales 
that characterize the fluid bridge are many orders of mag- 
nitude smaller than the macroscopic dimensions of the 
flow. Such a separation of length scales often leads to 
universal behavior as in the reverse process to coales- 
cence, drop break up, in which one drop breaks up into 
two or more droplets [ij |2] 13 [H |5] . In this case, the tran- 
sition proceeds as the radius of the liquid neck connect- 
ing the two pieces of fluid approaches zero. This causes 
the dynamics to approach a singularity, motivating the 
comparison of such fluid transitions to critical thermody- 
namic phase transitions [6, 7^ 8 . We also expect singular 
behavior during drop coalescence. However, recent dis- 
coveries have shown that not all fluid singularities obey 
universal dynamics [9, fTO] so that fluid transitions have a 
broader range of behaviors than the analogy with critical 
phenomena might suggest. 

For low- viscosity fluids such as water, the widening of 
the bridge during coalescence is opposed primarily by in- 
ert ial rather than viscous forces. Thus, throughout much 
of coalescence, the bridge radius r is much larger than the 
viscous length scale of water, Ij^ = ja'^ / ~ 14 nm. A 
straightforward scaling argument has been used to de- 
scribe the process pT] . 

The film of air between the drops ruptures to form a 
fluid bridge only when the drops are separated by a small 
distance. Soon after the bridge is formed, the length of 
the bridge, (i, is far smaller than its radius, r. Thus, d is 
the relevant length scale for calculating the pressure due 
to interfacial tension. Comparing interfacial tension to 




FIG. 1: Sequence of images showing the formation of a fluid 
bridge between two drops of aqueous NaCl solution at sat- 
uration, fluid density p = 1.1972 g/cm, kinematic viscosity 
ly = 1.662 cSt, and surface tension 7 = 82.55 dyne /cm [19). 
The drop radius A = 1 mm. The frames are separated by 69 
ps. 



inertia, we find 



^pd^ dt' 



(1) 



For a bridge spreading between two static, hemispheri- 
cal drops, d = I where A is the drop radius. The 
solution of this differential equation gives: 
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where to is the instant of coalescence and r = t — to. 

Simulations studying the coalescence of low-viscosity 
fluid drops in vacuum have confirmed this scaling law|T6l 
|T7|. Also, experiments using high speed imaging at up 
to 10^ frames per second have observed r cx r^/^ for 
r > 10 /isj2l [131 Ell- However, the speed and the ge- 
ometry of the transition limit the range of imaging stud- 
ies so that times earlier than 10 ps were impossible to 
access. Here we present an electrical method that al- 
lows us to study the fluid bridge between two coalescing 
drops at much shorter times: r ~ 10 ns. An electrical 
method was developed by Burton et al., where a small 
DC voltage was used to measure the resistance of a drop 
of mercury during break upjTS]. Our method adapts that 
technique to an AC voltage. This allows us to measure 
both the time-dependent resistance and capacitance of 
two coalescing drops of a conducting ionic solution such 
as aqueous NaCl. We find a new regime when r < lOps 
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FIG. 2: Experimental Setup. Two nozzles of radius A are se- 
cured in line with one another with tips separated by 2 A. A 
drop of aqueous NaCl solution is formed on the upper nozzle 
using a microliter syringe. The lower drop is slowly grown un- 
til the two drops coalesce. Electrodes are secured in the mea- 
surement cell and connected to a Wheatstone bridge circuit. 
The upper left branch consists of known circuit elements {Rt 
and Ct), while each lower branch is the input impedance of 
an NI PCI-5105 high-speed simultaneous sampling digitizer 
(Ro) in parallel with the capacitance of the coaxial cables 
(Co). The impedance of the cell can be separated into three 
contributions added in series: Zeiectrodes, Z fluid, and Zcr. 
Zeiectrode cau be modeled as a frequency-dependent capaci- 
tance in parallel with an equivalent resistance due to charge 
transfer, as shown in the dotted circle [18]. The impedance 
of the coalescing region, Zcr, can be modeled as a capaci- 
tor in parallel with a resistor, as shown in the dashed square. 
The conductivity of salt water at saturation is a = 0.225 {Q- 
cm)-^JJ, 



inconsistent with the predictions of the scaling argument 
outlined above. 

In our experiments, two nozzles of inner radius A were 
aligned vertically as shown in Fig. 2. We filled the noz- 
zles with saturated aqueous NaCl solution and formed 
a hemispherical drop at the tip of the upper nozzle by 
injecting a known quantity of fluid. We applied a fixed 
amplitude AC voltage of frequency / across electrodes se- 
cured opposite the nozzle tips, then grew the lower drop 
at a fixed rate until coalescence occurred. As shown in 
Fig. 2, we used a Wheatstone bridge to measure Zceii as 
a function of time during coalescence. We sampled two 
voltages simultaneously, Vr and Vg. With the exception 
of / = 10 MHz, which was sampled at the maximum 
rate of 60 MHz, the sampling was done at 10/. Using 
Labview (National Instruments), we measured the ratio 
of the amplitudes, IK-I/I^sl and the phase shift between 



them, Acj) as a function of time. This allows us to calcu- 
late the real and imaginary parts of Zceii- We determined 
to ±l/f using the shift in at contact. 

This cell impedance can be separated into three con- 
tributions added in series, as seen in Fig. 2: 
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We model the electrode impedance Zeiectrode using equiv- 
alent circuit elements. For a voltage across the cell 
\yceii\ ^ 50mF, the charge transfer from the electrodes 
into the solution can be modeled as an equivalent resis- 
tance. Additionally, the polarization of the fluid around 
the electrodes can be modeled as a frequency-dependent 
capacitor in parallel with this resistance [18 . The fluid 
between the electrodes and the nozzle tips contributes an 
impedance Z fluids and the "coalescing region" between 
the tips of the nozzles contributes an impedance Zcr- 
We bring the nozzle tips into contact to measure 
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Zceii - Z closed' We modcl Zcr as a resistor Rcr (repre- 
senting the resistance of the liquid in the shape formed 
by the two drops and the neck between them) in paral- 
lel with a capacitor Ccr (representing the capacitance 
between the drop surfaces). These quantities depend 
on the geometry of the coalescing region and are time- 
dependent. Our experiment measures Rcr and Ccr- 

Rcr versus r is shown in Fig. 3 for A = 1 mm. 
The filled symbols show twelve independent coalescence 
events, where three events are taken at each of four fre- 
quencies / ranging from 10 kHz to 10 MHz. The open 
symbols show this data binned logarithmically and av- 
eraged. We find that the data is described well by: 
Rcr = oiT-^ + /3r-i/2 ^ ^ ^^^^ g^^^^ y^^^ ^ig. 3 
shows this fit with a = (1.23 ± 0.3) • 10"^, = 0.8 ± 0.2 
and (5 = 23 ± 15. 

The inset to Fig. 3 shows Ccr^ the capacitance of 
the system during the 5 /is before coalescence occurs for 
three independent coalescence events taken at / = 10 
MHz. For / < 10 MHz, the signal before coalescence 
is comparable to the system noise. Cinit is constant 
within error, and the dashed line shows the average value, 
Ccr = 1.3 ± 0.14 pF. After coalescence, once the bridge 
is formed, the capacitance is poorly defined. 

We predict Rcr and Cinit ^ the capacitance of the co- 
alescing region just before coalescence occurs, using the 
scaling argument of Eq. 2. We separate Rcr into three 
resistances connected in series, Rcr — Rupper-^ Rbridge-^ 
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lower- Rupper {Rlower 

) is the resistance between the up- 
per (lower) nozzle tip and bridge, Rbridge is the bridge 
resistance. 



We model the drops such that R 
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Rhemi whcrc Rhemi is the resistance of a hemisphere trun- 
cated by a plane parallel to the flat surface of the hemi- 
sphere. The plane intersects the hemisphere with radius 
rtr, as shown in Fig. 4a. The resistance of this shape can 
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FIG. 3: Resistance during droplet coalescence. Rcr versus 
T = (t — to). A = 1 mm. The drops approach one another at a 
rate of 0.0004 A/ms. The closed symbols show 24 individual 
coalescence events, six obtained at each of four measurement 
frequencies. The open symbols are the average of the closed 
symbols, binned logarithmically. The error bars reflect the 
spread in these measurements as well as systematic error due 
to inaccuracies in the measurement of Zeiectrode and due to 
the choice of to, the instant of coalescence. The solid line 
shows Rcr = 1.23 • 10"^)t"^ + 0.79t"^/^ + 23.3. The inset 
shows Cinit versus r. Data is shown for three independent 
coalescence events at / = 10 MHz. The dashed line is the 
average value of the data shown, Cinit = 1.3 pF. 



be calculated numerically using the electrostatics calcu- 
lation package EStat (FieldCo). We vary rtr over several 
orders of magnitude, and find Rhemi = l/4crrtr, where a 
is the conductivity of the fluid. 

For two hemispherical drops just touching at the tips, 
d = 2r'^/A. Therefore, as Rbridge ^ d/airr'^^ we find that 
Rbridge = 2/(cr7rA), a Constant. 

The scaling argument summarized by Eqn. (2) shows 
that r oc r^/^, and calculating Rcr as a function of r, 
we find 
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For aqueous NaCl solution in air and drop radius A = 1 
mm, Rcr = 0.97r~^/^+28.3. This result is in qualitative 
agreement with the data for r > 10 fis. However, for 
r <C 10 /is, Rcr ^ . This is incompatible with the 
scaling argument. 

We can account for this discrepancy with a slight mod- 
ification of the geometry. Eq. 2 was derived assuming 
d (X r'^. However, if the drop tips are not quadratic, this 
is no longer true. For example, if prior to coalescence the 
drop tips are slightly flattened out to a radius rfiat^ as 
shown in Fig. 4b, then at early times, d is constant and 
only at late times would we see the d (x r^. Assuming we 



FIG. 4: Two geometries for coalescence, (a) Two hemispher- 
ical drops of radius A coalesce, forming a bridge of maximum 
radius rtr and height d = 2rtr/A. (b) Two drops coalesce with 
a flattened tip of radius rfiat- Here, for a bridge of radius r, 
while r < rfiat, d — constant. 



remain in the inviscid regime, Eqn. (1) becomes 



(5) 



In order to calculate Rcr in the altered geometry of 
Fig. 4b, rather than Rhemii we model the resistance 
of a hemisphere with a hole of radius r at the center 
of the flattened region radius Vfiat- We vary r, the ra- 
dius of the rapidly widening bridge, over several orders 
of magnitude. Solving numerically, we find that the re- 
sistance of this shape is Rfiat = l/(4rcr). Rbridge can 
be estimated as above, and for d constant, we find that 
Rbridge ^ d/cFTir'^. We are able to find Rcr versus time 
by combining this with Eq. 5: 
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This enables us to estimate the separation of the flattened 
tips, by comparing the predicted prefactor for the 
term to our measurements. We find d = 200 ± 100 nm. 

In our experiments, there is a negligible contribution 
to Rcr from the term. Examining Eqn. (6), we see 
a crossover from to behavior at r = 7 ± 5 ns for 
d — 200 ± 100 nm. Thus, we should see no contribution 
from the term at the earliest times accessible by our 
experiments. 

These conclusions are based on measurements of Rcr 
versus r. The inset to Fig. 3 shows our measurements 
of the capacitance Ccr- However, the capacitance of 
two nearly-touching spheres is only logarithmically de- 
pendent on their separation Even a small error in 
our measurement of Cinit leads to enormous uncertainty 
in the calculated separation. Moreover, since Ccr repre- 
sents the capacitance of the entire cell, the contributions 
from outside the "coalescing region" must be subtracted 
to obtain the relevant capacitance Qnit, which represents 
the capacitance of only the drop tips [21]. We measure 
Cinit = 0.41 ± 0.14 pF. This value is consistent with the 
model of the distorted drop tips but cannot exclude the 
original assumption of hemispherical drops at the instant 
of coalescence |21j. 
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In conclusion, we have observed an unexpected asymp- 
totic regime in the coalescence of two drops. This regime 
becomes visible at r < 10 /is, which is earlier than the 
shortest times accessible by previous imaging studies. 
Our electrical method allows us to study times three 
orders of magnitude earlier than this. Previous exper- 
iments, theory and simulations indicated that drops co- 
alesce while maintaining shapes described by quadratic 
minima. However, our data is inconsistent with that pic- 
ture and suggests that the coalescence occurs at the inter- 
face between two slightly flattened hemispherical drops. 

The scaling arguments and simulations discussed ear- 
lier do not account for such a flattening, and an un- 
derstanding of this phenomenon might illuminate other 
physical questions, such as the origin of the thin film 
rupture that initiates coalescence [22]. Also, industrial 
applications of fluid drops, such as coating procedures, 
Inkjet printing, and mixing in microfluidics could be af- 
fected by the way in which drops coalesce at the smallest 
scales and earliest times. 

Very near the instant of coalescence, the small-scale 
flows are decoupled from the large-scale flows. Because 
of this, we expect the analogy between topological tran- 
sitions in fluids and critical phase transitions would be 
most accurate at these earliest stages. However, we find 
that the scaling in the drop coalescence transition is de- 
pendent on the geometry and deformability of the drops. 
This has no analog in critical phase transitions. 
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